ON THE ABSOLUTELY CONTINUOUS SPECTRUM OF DIRAC 
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Abstract. We prove that the massless Dirac operator in with long-range 
potential has an a.c. spectrum which fills the whole real line. The Dirac 
operators with matrix-valued potentials are considered as well. 



1. Introduction 

In this paper, we consider Dirac operator for the massless particle in the external 
field generated by the long-range potential 



H 



-ia-sj + V 



(1) 



Here 





(To 



0-1 



0-2 



CT3 



Matrices aj are called the Pauli matrices. Denote the exterior of the unit ball in 
M.^ by 17. We use a notation E for the unit sphere in R^. Consider H in the Hilbert 
space [L^(r2)]^. Assume that the elements of a self-adjoint 4x4 matrix-function 
V{x) are uniformly bounded in O. We also assume that V has certain canonical 
form after the spherical change of variables. There are many meaningful potentials 
that satisfy these assumptions [20]. Consider the self-adjoint operator !K, generated 
by the boundary conditions f3{x) = fiix) = as a; S S. We introduce a matrix 

" 1 " 

10 

0-10 

-1 

The main result of the paper is the following statement. 

Theorem 13.11 Assume V{x) = v{x)(3, where v{x) is a real-valued, uniformly 
bounded, scalar function satisfying the following condition 



\x) 



+ 1 



< CO 



(2) 



Then, aac{'K) = M. 

This result can be regarded as a PDE (Partial Differential Equations) analog of 
the celebrated results by Szego on polynomials orthogonal on the unit circle with 
£^ Verblunsky parameters |17[ 119) . If one considers the power-decaying potentials 
v{x) : \v{x)\ < C{\x\ + 1)~°'^~'^, then condition Q is satisfied for any e > 0. In 
this case, we also obtain an asymptotics for the Green's function. 
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For the massless Dirac operator, theorem 13.11 solves Simon's conjecture JHI for 
Schrodinger operators. Under more conditions on V, the spectrum of 3{ is purely 
a.c. on R [2 2) . One can easily construct an example when conditions of the theorem 
13.11 are satisfied and the rich singular spectrum occurs. In the one-dimensional 
case, the first result on the presence of a.c. spectrum for slowly decaying potentials 
is due to M. Krein 0. See also 01 1121 E] for the modern development. The 
existence of wave operators for the one-dimensional Dirac operator with square 
summable potential was proved in [7j- An interesting paper |l(Jj discusses Szego- 
type inequalities for Schrodinger operators with short-range potentials. We will 
use some ideas from ,1(1. For discrete multidimensional Schrodinger operator with 
random slow decay, the existence of wave operators was proved by Bourgain P). 

The structure of the paper is as follows. In the second section, we consider one- 
dimensional Dirac systems with matrix-valued potentials. Then, in the third part, 
we deal with a multidimensional Dirac operator. 

The following notations will be used. \M\ — V M*M denotes the absolute value 
of matrix A/, symbol 1 will often stand for the identity matrix or operator. As usual, 
Cq°(R"'") denotes the space of infinitely smooth functions (or vector- functions) with 
the compact support inside (0,cxd). xk{x) denotes the characteristic function of 
the set K. {x, y) stands for the inner product of two vectors x and y in . 

2. One-dimensional Dirac operator with matrix-valued potential 

In this section, we study the one-dimensional case. Let us consider the Dirac 
operator in the following form 



D 



V, V = 



'b —a 
-a b 



(3) 



-d/dr 
d/dr 

where a(r),6(r) are m x m self-adjoint matrices with locally integrable entries. 
This form of a general Dirac operator is called a canonical form (^Jj PP-48-50). 
The elements of the Hilbert space are / = (/i,/2)* with /i(2) € [L'^{R+)]"'. The 
boundary condition /2(0) = defines the self-adjoint operator 2). We start with 
an elementary spectral theory of D. Consider solutions of the following equation 



D 



5- 



= A 



5- 



,$(0,A) = 1,*(0,A) = 0,A e 



For any /, consider the generalized Fourier transform 



^^(A)- / $*(r, A)/i (r)dr+ / ^'*(r, A)/2(r)dr 



There exists non-decreasing m x m matrix- function cr(A), A e R (spectral matrix- 
valued measure) such that ([Ml; P-106) 



(|/i(r)P + |/2(r)P)dr= / F*{X)da{X)F{X) 



(4) 



dajX) 
1 + A2 



< oo 



(5) 
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The resolvent kernel Rz {r, s) of operator T> has the following form 

oo 

1 



Rz{r,s) 



X-z 



$(r, A)da(A)$* (s, A) $(r, A)dcr(A)** (s, A) 
■^\r, A)da(A)$* (s, A) ^'(r, A)dCT(A)^'*(s, A) 



(6) 



where the integral is understood in the distributional sense. Notice that 

oo 



(7) 





where the integral converges due to (jSJ. Assume a{r) — b{r) = 0, if r > i?. Then, 
we can always find the Jost solution F(r, A): 



DF = XF, F{r, A) 



" F,{r,X) ■ 




—i 


. F2ir,X) 




I 



for r > A e C+ 



Notice that i^2(0,A) is an entire matrix- valued function. It is non-degenerate in 
C"*". Indeed, otherwise we would have the non-real eigenvalue for the self-adjoint 
operator D. The matrix 

" $(r,A) Fi(r,A) 
*(r,A) F2(r,A) 

is non-degenerate for r = 0, A G C+, therefore, it is non-degenerate for any r > 0. 
Consider the following matrix 

n -1 



Z{r,X) 



$(r,A) F^{r,X) 
vI'(r,A) i^2(r,A) 



Zii(r, A) Zi2(r, A) 
Z2i(r, A) Z22(r, A) 



Then, the resolvent of D can be written in the following form 



' hi ■ 




■ '^{r,X) ■ 








*(r,A) 


/ 



X)h2{s) + Zi2(s, X)hi{s)) ds+ 



F2(r A) / A)'^2(s) - Z22(s, X)hi{s)) ds, A e C+ 



Consequently, the resolvent kernel can be written in the form 



$(r,A)Zi2(s,A) -$(r,A)Zn(s,A) 
*(r,A)Zi2(s,A) -*(r,A)Zn(s,A) 



-Fi (r, A)Z22 (s. A) Fi (r, A)Z2i (s. A) 
-F2(r,A)Z22(s,A) F2(r,A)Z2i(s,A) 

In the free case {a — b ~ 0), $(r. A) = cos(Ar), ^'(r. A) = - 
— iexp(iAr), _F2(?', A) — exp(iAr), cr(A) = dX/i:, and 

i cos(Ar) exp(iAs) — cos(Ar) exp(iAs) 
—i sin(Ar) exp(iAs) sin(Ar) exp(iAs) 



Rl{r,s) 



i exp{iAr) cos(As) —i exp(iAr) sin(As) 
— exp(iAr) cos(As) exp(iAr) sin(As) 



, if r < s 

, if r > s 
sin(Ar), Fi{r,X) 

, if r < s 

, if r > s 



Notice that Z{0, A) 
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and 



Imi?A(0,0) = -[i?A(0,0)-i?^(0,0)] = 
Zi 



Ae C+ 



^ i_ \ f;-\o,x)f,*{o,x)-Fi{o,x)f^\o,x) o 

2i[ 

For A e C, let us introduce the foUowing functions: /i(r, A) = exp(— jAr)Fi (r, A), 
/2(r, A) — exp{—iXr)F2{r, A). Functions /i and /2 satisfy the foUowing equations 

/{ = -*A/i + a/i + A/2 - 6/2 
= -A/i - &/i - zA/2 - ah 

Therefore, a simple algebra yields 



(8) 



dr 

Using the asymptotics at infinity, we obtain 
A*(r,A)/2(r,A)-/2*(r,A)/i(r,A) 



-2ImA|/i+z/2 



(9) 



2 + 2ImA / \h{p,X)+if2{p,X)\'dp (10) 



Lemma 2.1. The following relations hold 

F^jr, X)F2{r, X) ~ Fj {r, X)Fi{r, X) 
i 

F*{0,X)F2{0,X)-F*{0,X)Fi{0,X) 



2,A e 



> 2,A e 



(11) 
(12) 



Proof. It suffices to use pU|l and definition of /i and /2. ■ 

Notice that the lemma is wrong if the Dirac operator is not in the canonical 
form. (|ll|l implies that ^2(0, A) is non-degenerate on M also. 

Corollary. The following relation holds: Imi?A(0, 0)(i4) ^3*^^(0, A)^^2"^(0, A), 
if A is real (see Therefore, using Q, we obtain 

<j'{X)^7r-'Fr\0,X)F,\0,X) (13) 

The factorization results of this sort are quite common in the scattering theory. 
Direct analogs hold for the so-called Krein systems (factorization via H-functions 
[3]), for one-dimensional Schrodinger operators with matrix- valued potentials (fac- 
torization via Jost functions Jj), for polynomials with matrix- valued coefficients 
(factorization via Szego function [SlEl^j). Notice that the Krein system with 
coefficient —A{r) yields operator —D 

Using the representation of resolvent by the solutions, one gets 

^F,{r,X)Fi\0,X) 



i?A(r,0) = 



-F2{r,X)F^H0,^) 



, for A e 



(14) 



So, we can find ^(0, A) by following a simple rule that turns out to be applicable 
to PDE 

-^2-1(0, A) 



lim R\{r, 0)e 



A e 
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Lemma 2.2. Consider a{r),b{r) g C^(M+). Then, for the function ^2(0, A), the 
following inequalities hold true 



Cy-' [\\a{r)r + \\h{r)r]dr 



\F2{Q,^V)\\ < exp 



The asymptotics 0/^2(0, A) as A e C+, |A| ^ 00 is 



i^2(0,A) = l-^ / \b{r) + za{r)\'dr + Oi\X\-' 



(15) 



(16) 



Proof. Consider a matrix 

-1 

1 

Upon multiplying by J, the equation DY = XY can be rewritten as 



J = 



Y' + QY^ -XJY, Q = -JV 



Take 



Ya = 



e'^'' -ie-'^'' 



and find the solution to (|17|l in the following form Y = YqS. 
Thus 

(-6 + ia)e~2*^'' 



S' ^^Y^-'QYoS^ 
Imposing condition 

on the 2m x m matrix 



(—6 — ia)e 



S'(oo, A) 



2i\r 

1 








Si 

S2 



we get F = FoS". System (UHJ can be rewritten in the following form 

00 

Si{r, A) = 1 + / e-^'^'{b{s) - ia{s))S2{s, X)ds 



(17) 



(18) 



(19) 



S'2(r,A)= / e'^'^'{b{s) + ia{s))Si{s,X)ds 



(20) 



These two equations are easy to iterate. We obtain the following integral equation 
for 



Si{r, A) = 1 + / e-^'^'ibis) - ia{s)) I e^'^\b{t) + ia(t))Si{t, X)dtds (21) 



2iXti 
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One needs to use Gronwall's inequality to obtain the first statement of the lemma. 
Iterating H21|l and integrating by parts, one proves H16(l . ■ 

Remark. Notice that the usual scattering coefficients A{X) and B{X) can be 
calculated by the formulas 

A{\) = l{F2{0,X) + iF,{0,X)) = 5i(0,A), 

1 2 (22) 
BiX) -(F2(0,A) -zFi(0,A)) = -182(0, X),Xe C 

If a, 6 e C^(R+), then the asymptotics of A{X) is 

00 

A{X) = 1 - -L / \b{r) + ia{r)\^dr + 0{\X\-^),X ^ 00, A e C+ (23) 
2iX J 


Using Hll|l . we get the conservation law 

|A(A)|2 ^1 + \B{X)\\ AeM (24) 

From (O, 

|A(A)|2 > 1 + |S(A)|2, AeC+ (25) 

Notice that if ^ £ C™, ^ 7^ 0, then the function In ||F2"^(0, A)f || is subharmonic 
in C+. Thus, the following inequality holds 



7rln||F2 {0,iym<y / ^^-^^ dX ^ - J ^^^—5 dX (26) 



2 

following theorem is straightforward 



Here we used asymptotics of (0, A) as |A| ^ cx) and its continuity in C+. The 



Theorem 2.1. Assume that | |a(r)| |, | |6(r)| | e L'^{R+), ^ e C", ||C|| = 1. T/ien, 
the following inequality holds 

oc 00 

/ i^i^^fj^^dA > / [\Hr)r + ||6(.)f ] rf. (27) 







Proof. Assume a{r),b{r) e Cj^(M+) first. Then, the estimate (|15|l from lemma 
12.21 implies 

\F,-\0,^y)\>e^p [||a(r)f + ||6(r)|p] dr j (28) 

So, H26|l gives the needed estimate. Now, consider any ||a('')||, ||&(?')|| G X'^(M+). 
We can find the sequence a„(r),6„(r) G Ci^(R+) such that 

{\\an - af + \\bn - b\f)dr ^ 







as n — > 00. Let cr„ denote the spectral measure for an,bn- Then, the second 
resolvent identity yields the weak convergence of measures cr,i to a. For each indi- 
vidual fT„, we have the estimate ()27ll . Then, using the standard argument involving 
semicontinuity of the entropy |H] , one gets . ■ 
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Remark. Constants in the inequality H27|l do not depend on to, the size of 
the matrix. Therefore, one can easily prove an analogous result for the Dirac 
operator with square summable operator-valued potential. As a simple corollary. 

Remark. The similar results can be proved for the Krein systems E| with 
matrix-valued or operator-valued coefhcients. Notice that the Krein system with 
coefficient —A{r) in a standard way generates the Dirac operator Dg = —T) with 
D introduced by ||2Jl. 

Assume that a(r), h(r) are integrable with finite support. Consider a fixed vector 
^ G C™. Let us find an element of the Hilbert space with generalized Fourier 
transform (A — iy)^^£^. This function is given by the formula 



(A-zy)-i$(r, A)da(A)e 



and is equal to i?ij,(r, 0)[^, 0]* due to ©. At the same time, (|14ll says that this 
vector is equal to 

F2{r,iy)F^\Q,iy)S, 
Therefore, the following arguments are valid. If the vector-function 

/i(r) = (/ii(r),/i2(r))*eC2'" (29) 

satisfies equation Dh = iyh and decays at infinity, then hi{r) ~ Fi{r, iy)'i], h2{r) = 
F2{'r,iy)ri, where 77 e C™. Take ^ = F2{0,iy)ri = ^.2(0). The generalized Fourier 
transform for the function h{r) is —(A — i?/)~^/i2(0). The spectral measure of h is 



a{X,h) 



d(a(A)fe2(0),fe2(0)) 
A2 + y^ 



Using the estimate H26f) . we have 

^ In [7r{<j'{X,h){X^+y^))] 
A^ + 



dX > 27rln||F2"^(0,i2/)/i2(0)|| = 27rln||77| 



or 



y I i^^f^dA>C + 2vrln|H 
X'' + 



(30) 



where the constant C depends on y. This inequality will play the crucial role later. 
Roughly speaking, it means the following. 

// we found at least one solution having the "free" asymptotics at some point in 
C^, then the entropy of the spectral measure can be controlled by the "amplitude" , 
i.e. \\t]\\. 
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The definition of the spectral measure gives an equality 

a'{X,h)dX = 



Therefore, the lower bound for \ \ri\\ yields the lower bound for 

oo 

ln-cr'(A» 



A2 + y2 



dX (31) 



Consider a(r),5(r) e Cg='(M+), 1 1^| | = 1. The functions In (0, A)C^OMn |(A-i(A)C, 01, 
ln||j4(A)|| are subharmonic in C+ and tend to as |A| oo, A € C+. Write the 
corresponding inequality at A = iy. Taking y — > cx), one gets the following estimates 
by comparing the coefficients against y~^. 

oo oo 

J ln(7rCT'(A)00dA = 2 J In 11^2-^(0, A)e|MA > 

— OO — OO 

oo oo 

>2 / \n\{F-\0,X)tO\dX>-C J [| |&(r)e + ^a(r)f | P] dr; 



oo 

\n\iA-\X)^,0\dX>~C [ [\\bir)^ + za{rm']dr; 



\n\\A-\X)\\dX >-C j [IkWlP + \\b{r)\\^]dr 

-oo 

Notice that estimates < 1, < 1 imply 

OO 

\n\{A-\X)tO\dX>-C I [\\b{r)^ + ta{rm'] dr 
^ , ^ (32) 

\n\\A-\X)\\dX >-C j [\\a{r)f + \\b{r)f]dr 

A 

for any A C K. Consider Dirac operator (Q. Fix some ^ £ C™, ||^|| = 1. Estimate 
suggests that condition (6(r) + ia{r))^ e i^(R+) guarantees aacCD) — R. 

Theorem 2.2. Let a{r) and b{r) have locally summable entries. Assume that there 
is a fixed vector ^ G C™ such that (6 + ia)^ G L'^{R+). Then, aad'D) = K. 

Proof. Without loss of generality, assume ^ — ei = {1,0, . . . ,0)*. Let us start 
with the compactly supported a and b. In this case, the Jost solution F(r, A) makes 
sense. Consider two functions /i(A) = A{X)F2^{0, A)ei and /2(A) = B{X)F2^{0, A)ei. 
From ll/i(A)f = \\F2\O, A)eif + ||/2(A)f for real A. Since /i + /a = e, (by 
1211), we have |(/i (A), ei)|2 = 11^2-1(0, A)ei||2+|(/2(A),ei)|2> 1(^2-1(0, A)ei,ei)|2 + 
|(/2(A),ei)|2. Consider 

(/i(A),ei) .... (/2(A), ei) 



(F2-i(0,A)ei,ei)"^' {F2\Q,X)e,,ei) 
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Then, 



a(A)+/3(A) 



(F2-i(0,A)ei,ei) 
and 

|a(A)|2 > 1 + |/3(A)|2,AgM 

Therefore, 

1(^2-1(0, A)ei,ei)| > -^,AeR 
2|a(A)| 

Function (i^2~"'^(0, A)ei, ei) is analytic in C+ and tends to 1 as |A| oo, A G C+. 
The numerator of a{X), function (/i(A), ei), does not have zeroes in C+. Indeed, in 
CT, we have ||/i(A)||2 > \\F^-\0, X)e4' + \\f2{X)\\' (by |23) and /i(A) + /2(A) = ei. 
Thus, |(/i(A),ei)|2 > 11^^2-1(0, A)ei||2 + |(/2(A),ei)|2 > \\F,-\0, X)^^ > 0. So, the 
function In |q;(A)| is superharmonic and nonnegative on M. For any A C M, we have 
the estimates 

J ln[7r(cr'(A)ei,ei)]dA > 2 J In ||F-i(0, A)ei ||dA > 2 J In |(F-i(0, A)ei, ei)|dA 



>-'2 J hi|a(A)|rfA-2|A| hi 2 

A 

Superharmonicity of ln|Q;(A)| imphes 



(33) 



7r\n\a{ty)\>y / ^^^dX (34) 



Consider a^jbe € C^(R+) that approximate a, b in L'^(R+) norm as e ^ 0. Denote 
the corresponding spectral measure by a^- Apply (|33|) and (|34|l . Take y — s- +oo in 
(1^ and use lemma 1^ and (f^ . We obtain 

oc 

ln|a,(A)|dA< I J \\{be{r) + iae{r))eifdr (35) 



Estimates (O and yield 

oo 

ln«(A)ei,ei)dA>Ci+C2 j \\{b,{r) + ia,{r))eifdr (36) 



Since dcTe converges weakly to du, we have 

CXD 

ln(a'(A)ei,ei)dA>Ci+C2 / \\{b{r) + ia{r))eif dr (37) 



Now, take arbitrary air) and b{r). Consider truncations a„(r) = a(''')x[o,„] ("r), 
bnii") = b{r)x[o^n]{'>')- Since the functions a„ and 6„ are compactly supported, 
estimate (pTzjl holds for the corresponding measure ct„(A). It is the general fact of 
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the spectral theory, that an converges to a weakly as n — > oo. Take n ^ oo to 
obtain 

oo 

J ln((T'(A)ei,ei)dA> C1 + C2 J || (6(r) + m(r))ei Ipdr (38) 

A 

That completes the proof. ■ 

Since the constants Ci and C2 in 1)38(1 are independent of m, the size of the 
matrices, the theorem is true for the operator-valued Dirac systems as well. Its 
statement is very strong, it proves a certain rigidity of the Dirac operator 
That makes theorem 12 . 21 applicable even to some PDE. 



3. Multidimensional Dirac operator 
In this section, we consider two operators 

H = -ia -XJ + V, Hs = -ia • V + 

The first one - on [L^(fi)] ^ with boundary conditions /s = /4 = on E. The second 

operator - on [L^(M'^)]''. Potential V{x) is always assumed to be symmetric 4x4 
matrix with the uniformly bounded entries. For iJ, it is given on fl, for Hg - on R'^. 
Thus, we have two self-adjoint operators !K and Jig- We also assume that after the 
"spherical change of variables"^, the matrix of V has the canonical form in a sense 
of section 2. The typical example of such a potential is V{x) — v{x)P, where v{x) 
is a scalar real- valued function (see 20 , p. 108 for other interactions and physical 
explanations). For simplicity, we will deal with this type of potentials only. 
For a, /3, we have 

akai + aiak = 2dki, auP + Pa^ = 0, = 1; fc, / = 1, 2, 3 

Notice that by letting / = — jcti, ./ = — icr2, K = — ids, we have relations 

/2 = J2 = = -1, 

IJ = - JI, IK = -KI, KJ = - Jif, IJ = K, JK = 1, KI = J 

which makes /, J, and K quaternions. The following algebraic relations are easy to 
verify. 

Lemma 3.1. If j E S, then 

(q!i7i -1- Q!272 + asTs)^ = 1 

(q!171 -I- Q!272 + Q!373 + 1)^ = 2(q!i7i + Q!272 + "373 + 1) 

(ai7i + 0272 + 0^373 + l)/3(ai7i + a272 + ^373 + 1) = 



Lemma 3.2. The resolvent kernel of the free Dirac operator has the following 
form p. 39) 

Gli., s) . ( .^^li^ + A^^^ + a) (39) 
V \x — sr \x — s\ J 47r x — s 



This rigidity is also present for polynomials orthogonal on the unit circle. 
'We will define this change of variables later. 
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Proof. Use the identities 



(J{^ + A)(J{r-A2)- 



-A 



(40) 



and the exphcit formula for the resolvent kernel of the free Laplacian. 
We will start with the following theorem. 



Theorem 3.1. Assume V{x) — v{x)f3, where v{x) is a real-valued, uniformly 
bounded, scalar function satisfying the following condition 

v'^{x) 



1x12 + 1 



< oo 



(41) 



Then, Uac.{'^) 



Proof. As in ^D], we can assume without loss of generality, that v{x) — in 
1 < |a;| < 2. Indeed, otherwise we subtract function f (a:)/3xi<|a;|<2(2;) from v{x)j3. 
The resolvent of JC" is an integral operator. One can obtain an expression for its 
kernel by using identity (|40|) with JC" instead of Jf^ . Therefore, the resolvent of "K is 
an integral operator too and the trace-class argument 13 can be applied. Now, we 
take any nontrivial infinitely smooth radially-symmetric function f{x) with support 
in {1 < |a:| < 2}. We will show that the spectral measure of the element 



f{x) = (/(r), 0,0,0)* 



(42) 



has an a.c. component which support fills R. The proof consists of two steps. In 
the first step, we represent operator "K in the canonical form ((SJ with unbounded 
operator-valued coefficients. Then, we apply the results of section 1 to obtain the 
needed estimates for the entropy of the spectral measure of /. 

Let us start with the suitable change of variables. Consider the Dirac operator 
!K. Let us write this operator in the spherical coordinates (see J20^, p. 126 or I21j . 

p. 16). The standard unitary operator f{x) e ^^(il) ^ F{r) = rf{ra),a S S maps 
scalar functions of three variables to a vector-valued function of one variable. For 
almost any r > 1, F{r) G L'^{'S). Following we consider the following auxiliary 
operators 

3 



-id/dxj 



p^ = -i[d/dr + r '^), Lj = Xj+iPj+2 



Xj+2Pj+l, S = CTO + E CTjLj 

where indices j are understood modulo 3, as usual. Notice that the operator pr is 
unitary to ^id/dr under U and s is independent of r. We want to get a matrix 
representation of !K convenient for us. To do that, we take a unitary, independent 
or r matrix U 

' <Jo 

—iffr 



u = 



Then, (im, p. 18) 

U*3{U = 



— ico 
iao 



s 
s 



-v{x)U*pU 



(43) 
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Operator s is selfadjoint in Denote its orthonormal cigcnfunctions by ip(n) 

where n is a multiindex. Therefore, in we can introduce an orthonormal 

basis spanned by the functions 



(ri) 







Then, any function f{x) G [L^(il)]'' can be represented as a sum /(x) ^ F{r) — 
^ (pj'^j (r)5'^jj + X] ('')^(fc) • The matrix of the unperturbed operator !Ko can 

be written as (^, p.l28 and (211, P-22) 



dr 








_d «(2) 

dr r 



_d_ 

dr 





_d "^(2) 

dr r 






where are eigenvalues of s. The boundary conditions at r = 1 are ip^^^{l) — 0. 

Now, we want to work with the new Hilbert space the space of vector- 
functions <i>(r) = {f^i-j{r), . . . , (/3^^(r), . . .)*, with the norm given by the formula 
below 



oo 



dr < 



So, the free multidimensional Dirac operator can be represented as an infinite or- 
thogonal sum of one-dimensional Dirac operators. Using elementary properties of 
the Pauli matrices, we get 

v{x)U*(3U = v{x)P 

Consider the multiplication by v(x) in [L^(S])]^ as a self-adjoint bounded operator 
with matrix —b{r), b(r) — {bij{r)},i, j = 1, . . . ,oo. Since the section 2 deals with 
an interval [0, oo) rather than [l,oo), we shift the argument r by one. We end up 
with the canonical representation (|3J), where a(r) is an unbounded operator in the 
diagonal form for each r > 0. Notice that the constant function (1,0)* G [L^(I])]^ 
is an eigenfunction of the operator s corresponding to the eigenvalue 1. For 5*^^ 

and ^("i), we choose the normed vectors coUinear to (1,0,0,0)* and (0,0,1,0)*, 
respectively. We have = 1 and the function /, given by H42|l . corresponds 
to $(r) — {{r + l)f{r + 1), 0, . . . , 0, . . .) in the vector-valued representation. From 
now on, we will deal with the representation (O of an operator CK. 

In the second step of the proof, we implement the main result of the second 
section, theorem 12.21 We can not do that directly, because we are dealing with an 
operator-valued coefficient a{r) unbounded for each r > 0. To avoid this difficulty, 
we consider operators D„ with bn = PnbPn, where Pn denotes the orthogonal 
projection onto the linear combination of the first n functions V'(fe) ■ The matrix of 
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bn{r) is 



bn{r) 



bii{r) 

bni{r) 




binir) 

bnn{r) 




,r > 



Each of 2)„ can be written as an orthogonal sum of the Dirac operator !)„ with 
matrix- valued potential Vn of size nxn and an infinite number of the one-dimensional 
Dirac operators with scalar potentials. Denote by iL„ the subspace of H on which 
T>n acts. Matrix Vn has the following form 



with 



a„{r) 



-bn{r) -a„(r) 
-an{r) bn{r) 



,r > 



^ 

• • • 







r+1 



Notice that the function <I>(r) lies in for any n. Denote the spectral matrix- 
valued measure of I)„ by ct„(A). Since b{r) = on [0,1], the spectral measure 
dfin{X) of $(r) is equal to |p(A)p((icr„(A)ei, ei), where 



p{\) = y (r + l)/(r + l) 





cos(Ar) 



sin(Ar) 



Here, 



, sin(Ar) rcos(Ar) , , 

cos(Ar) + — \; -sin(Ar) 1 + 



dr 



1 



A ' A(r + 1) —^--^ y ■ A2(r + 1)^ 
are generalized eigenfunctions for the Dirac operator Q with potential 

1 



V = 





1 

r + 1 



r+1 




From the proof of theorem 12. 21 festimate pTfl ). we get 

oo 

ln«(A)ei,ei)dA>Ci + C2 J || (6„(r) + ia„(r))ei H^dr 

A 

for any finite interval A g R. Notice that 

oo 

\\a„{r)ei\\^dr ^ J {r + l)-^dr 





(44) 



and 



|&n(r)ei|r < ||6(r)ei||^ = J v'i{r + l)T)dT 
res 
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So, 

OO oo 



\\bn{r)eifdr< J J v^^r + l)T)dTdr < C J j^^dx 
res n 

and we clearly have 

j ln/Lt;(A)dA > C > -oo (45) 
A 

with some constant C independent of n. Notice that bn b strongly. So, dfin{X) 
converges weakly to d/x, the spectral measure of / with respect to the initial operator 
Ji. The semicontinuity of the entropy |S] implies 



J \n^l'{X)d\ >C>-oo 



A 

■ 

The reduction of !H to a one-dimensional system with the operator- valued poten- 
tial is, probably, not necessary. One could have introduced the radiative operator 
and worked with this operator directly avoiding an approximation by Dirac opera- 
tors with matrix- valued potentials (see In the meantime, this reduction is not 
too difficult and it shows how very general facts for the matrix-valued orthogonal 
systems are applied to different PDE's. 

We do not consider the question of existence of wave operators. It might be that 
the problem can be solved by using an approach of W\. 

In the next theorem, we establish an asymptotics of the Green's function for the 
operator 'Kg. For JC, that means existence of the function h (see (I29II ') satisfying 
homogeneous equation Hh = Xh, X G C+ with the well-controlled "amplitude" . Due 
to H30|l and (P^, we have a certain estimate on the entropy of the corresponding 
spectral measure. That, in particular, gives another proof of aacC^) = R for the 
case of a power-decay. 

Theorem 3.2. Assume v{x) is given on MP and satisfies an estimate \v{x)\ < 
Cv{\x\ + 1)~*''^~'^, with fixed e > and sufficiently small. Then, the resolvent 
kernel Gx{x,0) of operator Hg at a point X — i has the following representation 



G.(x,0) = ^ 
47r a; 



i^+i)?i(x) + (|a;| + l)~"-^T2(x) 



(46) 



where 



|J'2(:^)|<C (47) 

uniformly in , || CPi — 1|| < (5, 5 ^ as C„ 0. Positive constant C depends on 
Cy and e only. 

We need the smallness of to guarantee convergence of a certain series. In 
general situation, one can take the spectral parameter A sufficiently far from the 
spectrum. Or, we can take X — i, and let potential be zero on the large ball centered 
at origin. That would make constant as small as we want^ but would not change 
the scattering picture. 

Let us prove some auxiliary lemmas first. 



With respect to a smaller e. 
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Lemma 3.3. The following estimate holds (1 < p < 2|a;|/3) 

J e-^^'-y^-^y^dTy < Cpe-\''\ 



(48) 



\v\=p 



\y\=p 



^-\x-y\-\y\ -2 



C{x,y)dTy < Ce 



\y\=p 



where (^{x,y) is an angle between x and y. 

Proof. Without loss of generality, assume that x — (0,0, Introducing the 
spherical coordinates j/i = p cos 9 cos Lp , y2 — p cos sin 93, 1/3 = ps\n9, we get 

p^ j dip j d6 cos6exp{~p - + p'^ - 2|a;|psin6') 



-7r/2 
7r/2 



< Cp'e 



d6'cos6'exp [-c|x|/9(|x| - p) ^l-sin^)] < Ce 



\x\ \x\ - P 



Estimate H48|) is now straightforward. The other statement of the lemma can be 
proved similarly. ■ 

Take any two vectors x,y e R^. The following inequality is obvious 



y_ ^ 

\y\ \x\ 



< C sin C 



(49) 



Together with ({x,y), consider xi^jU) ~ angle between x — y and a;. From the 
sine-theorem, sinx = \y — x\~^\y\ sin^. Consequently, 



X — y X 



\x-y\ \x\ 
By the triangle inequality. 



< Csinx = C|y — a;| ^|y|sin(^ 



For \y\ < 2\x\/3, 



y 


X — 


y 


< 


y 




X 


+ 


X 


x-y 


\y\ 


\x- 


y\ 




\y\ 




~\ 




R " 


\x - y\ 






y 




X — 


y 


< CsinC 








\y 




X — 


y\ 











(50) 



(51) 



(52) 



Let |a;| > 1 and T = {y : \y\ > 2|x|/3, \x - y\ > 2\x\/S}. 
Lemma 3.4. The following estimate holds 

exp(-|a; - y\ ~ \y\)dy < C exp(-7|x|) 



(53) 



with 7 > 1. 
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Proof. Indeed, in T, 

\x-y\ + \y\ > |2;|/5 + l6|x|/15 

Taking 7 = 16/15, we obtain the statement of the lemma. H 

In the following three lemmas, we will be estimating certain integrals over the 
M.^. Lemma 13.41 shows that the contribution coming from the integration over T is 
small and can be neglected. 

Lemma 3.5. The following bound is true 

r g-\^-y\ g-\y\ g-l^l 

L \x-y\' |y|1.5+^ + l^^ < ^|x|l-5 + l 



Proof. By lemma IX^ 



\x\ 

'v\ g-\y\ r g-k-al g-|y| 



|y- 



< ^ / < c, 







Lemma 3.6. The following relation is true 

/n- . (t. - lA \ p.~'\y\ 



f e-\--y\ f a-{x-y) ^ \ 
Jw \x-y\\ \x-y\ ) 



-dy 



|yP+e + 1 
+ l)^i(x) + (|a;| + l)-"-V2(x) 



\x\ + l 

where ipi(2)ix) are matrix-functions uniformly bounded in R"^. 
Proof. 

e-l""^l f a-jx-y) \ e-\y\ 
.\x-y\\ \x-y\ ^ ) I2/P+- + I ^ 
e^l^^^l / a ■ (x — y) a ■ x\ e"'^' 

k-yl V \x~y\ \y\^+' + i ^ 

a-x \ f e^l^^^^l e^l^l 

'1 / ^-^^^T^'^y = h+h 



\x\ J 7r3 \x-y\ \y\^+' + l 
Following the proof of lemma we get 



r g~\x-y\ g-|y| 

L \x-y\ |#- + i'^^<^RTT ^^^^-^P') 
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For Ii, we have 



a ■ (x — y) a ■ X 



x-y\<2\x\/3 Ix-yl V \x-y\ \x\ J + 1 



rdy 



-\y\ 



a ■ y a ■ X 



\x-y\ 



'\v\<2\x\/3 \y\ V \y\ \^\ J \x-y\^+^ + i 

By (gHJ and lemma lO 

\x\ 



dy 



|y|<2b|/3 



\y\ 



-dy < 





We now estimate integral in 1\ over \y\ < 2|a;|/3. By H50|l and lemma 1531 

g-k^al / Q, . (^x — y) a ■ x\ e^l^'l 



5+e 



< 



y\<2\x\/3 



Lemma 3.7. If v{x) is a real-valued scalar function satisfying an estimate 
\v{x)\ < {\x\ + l)^*''^^'^, then the following representation holds 

\x-v\\ \x-y\ ) \ \y\ ) \y\ 



\x\ + l 



+ lUi(a:) + (|x| + l)-"-V2(x) 



where ifn2){x) are, again, matrix-functions uniformly bounded in M"^. 
Proof. From the lemma ISTTl we infer 



Therefore, 



\x-v\ J \ \x-y\ 
g-\^-y\ ^ 01 . {x — y) 



,-\x- 



V - v\ 



x-y\ \ v-y\ ' \y\ \y\ 

a ■ (X — y} ^ \ „ , ^ I a ■ y a ■ {x — y)\ e~^y^ 



\v\ 



■ (x — y) \ „ , ^ f a ■ 



\y\ \x-y\ J \y\ 

J1 + J2 



where 



, a - X 

Ji - ( + 1 



g-\^-y\ 
\x - y\ 



■Pv{y) 



a ■ y a ■ (x — y)\ e 



\y\ \x-y\ J \y\ 



-dy 



and 



x~v\ \ \x-y\ 



\y\ \x-y\ J \y\ 
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Consider Ji. Using (|5^ and lemma rOl we obtain 



\x-y\ 



\y\<2\x\/3 



\x - y\ 



Hy)\ 



< c- 



-\x\ 



\x\ 



„0,5 



\y\ \x-y\ 



\v\ 



dy 



x\ + l J + 1 





dp < C- 



X +1 



\y\ 

(sharp!) 



\x-y\ 



y-x\<2\x\/3 



\x-y\ 



Hy)\ 



y x-y 



yl<2|a;|/3 



\x - y\ 



\v{x - y)\ 



\y\ \x-y\ 
y x~y 



-\v\ 



\y\ 



-dy 



\v\ 



< C- 



|1.5 + £ 



IJ p + l 



\y\ \x-y\ 

0.5 

^ dp<C 



\y\ 

-\x\ 



dy 



Using H5Q(I . (|52(l . and lemma rOl we get the following inequalities to estimate 



\J2\ 



\J2\< 



\x-y\ 



|y|<2|:r|/3 \x - y\ 



X — y X 



\x-y\ \x\ 



Hy)\ 



y x-y 



\y\ \x-y\ 



,-\y\ 



-dy 



< 



For the region jx — y| < 2|a;|/3, we have 



/ 

J \y-x 

L 



y-x\<2\x\/3 



\y\<2\x\/?, 



\x - y\ 
\x - y\ 



X — y X 



Hy)\ 



\x-y\ \x\ 

v{x - y)\ 



y x-y 



y_ x_ 

\y\ \x 



\y\ \x-y\ 
y x-y 



< 



Ce 



\y\ \x-y\ 

\x\ 



f,-\y\ 

Ha 

\y\ 



dy 
dy 



dp Ce-I^l 
< 



2;|1.5 + £ J |2;|1.5 ^ I 



Proof of the theorem l3.2l Let us iterate the second resolvent identity to get the 
needed estimate for Gi{x, 0) 

G.(x,0) = G°(a;,0) - / G^,{x, s)^3v{s)G^{s,Q)ds 



Using lemmas [3.51 - 13.71 and explicit formula for G^{x, s), we see that the n— th term 
in the corresponding series has the following form 

p-|a^l 



X +1 



a ■ X 



+ l]^[-\x) + {\x\+l)-'^-'^i-\x) 



(54) 



with \(p^^ll){x)\ < [C(£)a]". That can be easily proved by the induction. Indeed, 
for the first term, we have the representation H39|l . Assume that H54|l holds for the 
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n-th term T„. Then, 

Tn+i{x) ^-J G^{x, s)/3v{s)T^{s)ds = -{h+h+h + h) 

R3 

where 



J_ /■ g ■ (x -~ s) / exp(-|x-s|) \ ^ . X / _expH£|)_^^ („)^ 
Airj \x-sY \ \x-s\ )" + 



We apply lemma rTTI to /i, lemma rTHl to I2, and lemma PT^ to I3 and /4. If Cv is 
small enough, we will get uniform convergence of the series and the estimate (I47II . 
Clearly, || Ti - 1|| < (5 with 5 ^ if a ^ 0. ■ 

Remark. It is clear that the asymptotics of the Green's function Gi{x, s), \s\ < 1 
can be obtained similarly: it is close to Gl{x, s) as |x| — s- 00. 

Fix any £ > 0. Consider the Dirac operator IK with \v{x)\ < C„(|x| + 1)^°-^^^, 
where is sufficiently small, and operator 'Kg on M.^ with = on |x| < 1. 

Then, for any nonzero function f{x) — (/i, /2, /a, JaY G [L^(M'^)] with support in 
the unit ball, the function 



h{x) = J Gi{x, s)f{s)ds 

satisfies equation Hh = ih {h does not satisfy boundary condition on E unless h = 
in fl). Therefore, by the theorem 13 . 21 and by the remark above, we can control the 
"amplitude" of h. For |a;| large enough, 

e\-\\x\h{x)^^(^ + l] [ eM^,s)f{s)ds + f{x) 



Att \ \x\ J J \x 

and 

ll/»ll<<^ll/lli 

with 5 ^ as Ci, — > 0. For the integrals, we have 

X 



exp( — ,s)|/j(s)|ds > e ^ / |/j(s)|ds, (j = 1, . . . , 4) 



a; 



The "amplitude" is calculated as 

A{\x\-^x)^ lim \x\exp{\x\~l){h3{x),h4{x)y 

\x\ — >oc 

Take /i = /2 = and fs, - arbitrary non- negative. Then, A is well-defined 
as an element of [L^(E)]^ for Vn — v{x)x\x\<n{x)- Denote this sequence by An- 
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Now, theorem 13 . 21 savs that An does not go to zero. Using the reduction to a one- 
dimensional system and H3Q|I . one can show that the spectral measure of an element 
/ has an a.c. component which fills R. 

It is likely, that more careful estimates might prove some analog of theorem 13. 21 
for V satisfying condition One of the main ideas of the paper is as follows. Take 
the spectral parameter as far from the spectrum as we wish. If we manage to find 
at least one solution to homogeneous problem that satisfies certain asymptotics at 
infinity with nonzero "amplitude" , then an a.c. spectrum fills the whole line. The 
main advantage is that we need to do some rather delicate analysis of, say. Green's 
function far from the spectrum, not on the real line. That is an interesting problem 
to study the Green's function for operator |^ when the potential is random with 
slow decay or with no decay at all (the Anderson model). The approach of theorem 
13. II does not work in this case. In the meantime, analysis done in 3^ together with 
estimates used in the proof of theorem 13 . 21 seem to be relevant for the problem. 
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